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Context. Dense and cold clouds seem te populate the broad line region surrounding the central black hole in AGNs. 
These clouds could interact with the AGN jet base and this could have observational consequences. 
Aims. We want to study the gamma-ray emission produced by these jet-cloud interactions, and explore under which 
conditions this radiation would be detectable. 

Methods. We investigate the hydrodynamical properties of jet-cloud interactions and the resulting shocks, and develop 
a model to compute the spectral energy distribution of the emission generated by the particles accelerated in these 
shocks. We discuss our model in the context of radio-loud AGNs, with applications to two representative cases, the 
low-luminous Centaurus A, and the powerful 3C 273. 

Results. Some fraction of the jet power can be channelled to gamma-rays, which would be likely dominated by syn- 
chrotron self-Compton radiation, and show typical variability timescales similar to the cloud lifetime within the jet, 
which is longer than several hours. Many clouds can interact with the jet simultaneously leading to fluxes significantly 
higher than in one interaction, but then variability will be smoothed out. 

Conclusions. Jet-cloud interactions may produce detectable gamma-rays in non-blazar AGNs, of transient nature in 
nearby low-luminous sources like Cen A, and steady in the case of powerful objects of FR II type. 
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1. Introduction 

Active galactic nuclei (AGN) consist of an accreting super- 
massive black hole (SMBH) in the center of a galaxy and 
sometimes present powerful radio emitting jets (Begelman 
et al. 1984). Radio-loud AGNs have continuum emission 
along the whole electromagnetic spectrum, from radio to 
gamma rays (e.g. Boettcher 2007). This radiation basically 
comes from the accretion disc and bipolar relativistic jets 
originated close to the central SMBH. Radiation of accre- 
tion origin can be produced by the thermal plasma of ei- 
ther an optically-thick geometrically-thin disc under effi- 
cient cooling (Shakura & Sunyaev 1973), or an optically- 
thin geometrically-thick corona (e.g. Liang & Thompson 
1979). The emission from the jets is non-thermal and gen- 
erated by a population of relativistic particles likely ac- 
celerated in strong shocks, although other mechanisms are 
also possible (Rieger et al. 2007). This non-thermal emis- 
sion is thought to be produced through synchrotron and 
inverse Compton (IC) processes (e.g. Ghisellini et al. 1985), 
although hadronic models have been also considered to ex- 
plain gamma-ray detections (e.g. Mannheim 1993, Miicke 
& Protheroe 2001, Aharonian 2002). 

In addition to continuum radiation, AGNs also present 
optical and ultra-violet lines. Some of these lines are broad. 
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emitted by gas moving with velocities Vg > 1000 km s~^ 
and located in a small region close to the SMBH, the so- 
called broad line region (BLR) . The structure of this region 
is not well known but some models assume that the mate- 
rial in the BLR could be formed by dense clouds confined 
by a hot (T - 10^ K) external medium (Krolik et al. 1981) 
or by magnetic fields (Rees 1987). These clouds would be 
ionized by photons from the accretion disc producing the 
observed emission lines, which are broad because of the 
cloud motion within the SMBH potential well. An alterna- 
tive model proposes that the broad lines are produced in 
the chromosphere of evolved stars (Penston 1988) present 
in the nuclear region of AGNs. 

The presence of material surrounding the base of the 
jets in AGNs makes jet-medium interactions likely. For in- 
stance, the interaction of BLR clouds with a jet in AGNs 
was already suggested by Blandford & Konigl (1979) as 
a mechanism for knot formation in the radio galaxy M87. 
Also, the gamma-ray production due to the interaction of 
a cloud from the BLR with a proton beam or a massive 
star with a jet were studied in the context of AGNs by Dar 
& Laor (1997) and Bednarek & Protheroe (1997), respec- 
tively. 

In this work, we study the interaction of BLR clouds 
with the innermost jet in an AGN and its observable con- 
sequences at high energies. The approach adopted is similar 
to that followed in Araudo et al. (2009) for high-mass mi- 
croquasars (for a general comparison between these sources 
and AGNs see Bosch- Ramon 2008), where the interaction 
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of stellar wind clumps of the companion star with the mi- 
croquasar jet was studied. Under magnetic fields below 
equipartition with the jet kinetic energy (i.e. the jet should 
be matter dominated), cloud penetration will lead to the 
formation of a relativistic bow shock in the jet and a slow 
shock inside the cloud. Electrons and protons can be ef- 
ficiently accelerated in the bow shock and produce non- 
thermal emission, in situ via synchrotron and synchrotron 
self-Compton (SSC) mechanism, and in the cloud through 
proton-proton (pp) collisions. For magnetic fields well below 
equipartition, the SSC component becomes the dominant 
electron cooling channel, that leads to significant gamma- 
ray production. Since the bow shock downstream is almost 
at rest in the laboratory reference frame (RF), this emission 
will not be significantly boosted. The resulting spectrum 
and the achieved luminosities in one jet-cloud interaction 
depend strongly on the magnetic field, the location of the 
interaction region, the cloud size, and the jet luminosity. 
However, many clouds could be inside the jet simultane- 
ously, and then the BLR global properties, like size and 
total number of clouds, would also be relevant. Depending 
on whether one cloud or many of them penetrate into the 
jet, the lightcurve will be fiare-like or rather steady, respec- 
tively. 

In order to explore the radiative outcomes of jet-cloud 
interactions in AGNs, we apply our model to both Faranoff- 
Riley galaxies I (FR I) and II (FR II). In particular, we con- 
sider Centaurus A (Cen A) and 3C 273, the nearest FR I 
and a close and very bright fiat spectrum radio quasar (with 
FR II as parent population), as illustrative cases. Although 
in FR I the BLR is not well-detected, clouds with similar 
characteristics to those found in FR II galaxies may sur- 
round the SMBH (Wang et al. 1986, Risaliti 2009). Cen A 
has been detected at high- (HE) (Hartman et al. 1999; 
Abdo et al. 2010) and very high-energy (VHE) gamma rays 
(Aharonian et al. 2009), whereas 3C 273 has been detected 
so far only at HE gamma rays (Hartman et al. 1999; Abdo 
et al. 2010). We have computed the contribution of jet-cloud 
interactions to the gamma-ray emission in these sources, 
and estimated the gamma-ray luminosity in a wide range of 
cases. We find that gamma rays from jet-cloud interactions 
could be detectable by present and future instrumentation 
in nearby low-luminous AGNs at HE and VHE, and for 
powerful and nearby quasars only at HE, since the VHE 
radiation is absorbed by the dense nuclear photon fields. In 
the case of sources showing boosted gamma rays (blazars), 
the isotropic radiation from jet-cloud interactions will be 
masked by the jet beamed emission, which will not be the 
case in non-blazar sources. 

The paper is organized as follows: in Sect. [21 the dy- 
namics of jet-cloud interactions is described; in Sect. [31 a 
model for particle acceleration and emission is presented for 
one interaction, whereas in Sect. [3] the case of many clouds 
interacting with the jet is considered; in Sects. [5l and [51 the 
model is applied to FR I and FR II galaxies, focusing on 
the sources Cen A and 3C 273; finally, in Sect. [71 the results 
of this work are summarized and discussed. We adopt cgs 
units through-out the paper. 

2. The jet-cloud interaction 

Under certain combinations of the jet ram pressure, and the 
cloud size and density, cloud-jet penetration is expected to 
occur. The details of the penetration process itself are com- 





Fig. 1. Sketch, not to scale, of an AGN at the spatial scales 
of the BLR region. In the top part of the figure, the inter- 
action between a cloud and the jet is also shown. 



Table 1. 

jets. 



Values assumed in this work for BLR clouds and 



Description 


Value 


Gloud size 


Re = 10" cm 


Gloud density 


He = 10" cm"^ 


Cloud velocity 


Vc = 10^ cm s-i 


Cloud temperature 


Tc = 2 X 10'' K 


Jet Lorentz factor 


r = 10 


Jet half-opening angle 


(^«6° 



plex. Here we do not treat them in detail, but just assume 
that penetration occurs if certain conditions are fulfilled. 
For low magnetic fields, a cloud inside the jet may repre- 
sent a hydrodynamic situation in which a supersonic fiow 
interacts with a body of approximately spherical shape at 
rest. The cloud, as long as it has not been accelerated by the 
jet ram pressure up to the jet speed (fj), produces a pertur- 
bation in the jet medium in the form of a steady bow shock 
roughly at rest in the laboratoty RF, with a velocity with 
respect to the jet RF approximately equal to Vj. Since the 
cloud is not rigid, a wave propagates also through it. Since 
the cloud temperature is much lower, and the density much 
higher than in the jet, this wave will be still supersonic 
but much slower than the bow shock. The jet pressure ex- 
erts a force on the cloud leading to cloud acceleration along 
the axis, hydrodynamical instabilites and, eventually, cloud 
fragmentation. In the following, the jet-cloud interaction is 
described. Further discussion, and a proper account of the 
literature, can be found in Araudo et al. (2009). Sketchs 
of the jet-cloud interaction and the jet-BLR scenario are 
shown in Fig. [I] 

We adopt here clouds with typical density ric — 
10^° cm-3 and size i?c = 10^^ cm (Risaliti 2009). The veloc- 
ity of the cloud is taken Vc — 10^ cm s~^ (Peterson 2006). 
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The jet Lorentz faetor is fixed to F = 10, implying Vj ~ c, 
with a half-opening angle </) « 6°, i.e. the jet radius/height 
relation fixed to Rj = tan((/)) z = 0.1 z. All these parame- 
ters are summarized in Table [T] and will not change along 
the paper; from them, the jet density Uj in the laboratoty 
RF can be estimated: 



(F — 1) nip c'^cTj 



8 X 10* 



1044 erg s- 



10 



16, 



(1) 



where CTj = ttR? and Lj is the kinetic power of the matter 
dominated jet. 

The jet ram pressure should not destroy the cloud before 
it has fully entered into the jet. This means that the time 
required by the cloud to penetrate into the jet should be: 



ic 



2i?c 



2 X 10* 



Re 



10" cm/ VlO^cms-i 



(2) 



should be shorter than the cloud lifetime inside the jet. To 
estimate this cloud lifetime, let us compute first the time 
required by the shock in the cloud to cross it (tcs)- The 
velocity of this shock, Wcs, can be derived making equal the 
jet and the cloud shock ram pressures: (F — 1) nj nip c^ = 
ric TOp w^g, valid as long as Vcs '^ c . Then: 



X-^/^c^3xl0«f ,!"' , ) 
^ V10i"cm-3/ 

V 1016 cm/ 



10*4 erg s-i 



-1/2 



cms 



(3) 



where x is the cloud to jet density ratio, nc/nj{T — 1). This 
yields a clould shocking time: 



^ ^ 7 X 10* 



Re 



1013 cm/ VlOiOcm-3 



1/2 



UOiScnJ V 1044 erg s-iy 



-1/2 



(4) 



Therefore, for a penetration time (ic) at least as short as 
~ tcs, the cloud will remain being an effective obstacle for 
the jet flow. Setting ic ~ ics, we obtain a minimum value 
for X and hence for z. 

Hydrodynamical instabilities produced by the interac- 
tion with the jet material will affect the cloud. First of all, 
the jet exerts a force on the cloud through the contact dis- 
continuity. The acceleration applied to the cloud can be 
estimated from the jet ram pressure Pj, and the cloud sec- 
tion fjc ~ TT i?^ and mass Mc ~ (4/3) ttR^ nj uip: 



9 






4xi?c 



3 Vc 



(5) 



Given the acceleration exerted by the jet in the cloud, 
Ray leigh- Taylor (RT) instabilities will develop in the cloud, 
at the jet contact discontinuity, with timescale: 



(6) 



where the instability length I is the spatial scale of the per- 
turbation. For perturbations of the size of the cloud: / ~ Re, 





which are those associated to cloud significant disruption, 
one gets iRT ^ ^cs- 

In addition to RT instabilities, Kelvin-Helmholtz (KH) 
instabilities also grow in the cloud walls in contact with 
shocked jet material that surrounds the cloud. Accounting 
for the high relative velocity, Wroi ^ Vj, one obtains: 



(7) 



where ^i-ci ^ c^/x'- For ^ ~ ^c, we obtain ^kh ^ tcs- In the 
previous estimates of ipiT and ixH we have not taken into 
account the effect of the magnetic field (e.g. Blake 1972), 
since we assume that it is dynamically negligible. We note 
that, given g, the time to accelerate the cloud up to the 
shock velocity Wcs is '^ ics- However, the timescale to accel- 
erate the cloud up to Vj is ^ ics provided that Vj ^ Wcs, 
and before, the cloud will likely fragment. 

Finally, there are two additional timescales also rele- 
vant for our study, the bow-shock formation time, tbs, and 
the time required by the cloud to cross the jet, tj. The 
timescale ibs can be roughly estimated assuming that the 
shock downstream has a cylindrical shape with one of the 
bases being the bow shock, relativistic shock jump condi- 
tions, equal particle injection and escape rates, and a escape 
velocity similar to the sound speed ~ c/v3 (for a relativis- 
tic plasma). This yields a shock-cloud separation distance 
of Z ^ 0.3 i?,c, which implies: 



ibs 



Re 



- = 10^ , 

c V1013 



(8) 



Since in general tbs ^ ^cs, we can assume that the bow 
shock is in the steady regime. The jet crossing time ij can 
be characterized by: 



U 



2i?j 



2 X lO"" 



V 1016 cm/ V109cms-i/ 



(9) 



Note that if the cloud lifetime is < ij , the number of clouds 
inside the jet will be smaller than expected just from the 
BLR properties. 

In order to summarize the discussion of the dynamics 
of the jet-clump interaction, we plot in Fig. [2] the ics (for 
different Lj), tj, ic, and tbs as a function of z. As shown in 
the figure, for some values of z and Lj the cloud could be 
destroyed by the jet before full penetration, i.e. tcs < tc- 
This is a constraint to determine the height z of the jet at 
which the cloud can penetrate into it. Note also that, in 
general, tbs is much shorter than any other timescale. 



2.1. The interaction lieiglit 

The cloud can fully penetrate into the jet if the cloud life- 
time after jet impact is longer than the penetration time 
(the weaker condition that jet lateral pressure is < ric nip v'^ 
is then automatically satisfied). This determines the min- 
imum interaction height, Zjnt, to avoid cloud disruption 
before full penetration. Also, this interaction cannot oc- 
cur below the jet formation region, zq '^ 100it!g « 1.5 x 
lOi5(A'fbh/lO*A'f0) cm (Junor et al. 1999). For BLR-jet 
interaction and cloud penetration to occur, the size of the 
BLR should be i?bir > -Zo and Zint. 

The lifetime of the cloud depends on the fragmenta- 
tion time, which is strongly linked to, but longer than, tcs. 
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Fig. 2. The jet crossing (blue dotted line), cloud pene- 
tration (red dot-dashed line), bow-shock formation (violet 
dashed line) and cloud shocking (green solid lines) times 
are plotted; all of them have been calculated using the val- 
ues given in Table [T] The time tcs is plotted for Lj — 10"''*, 
lO'"^ and 10*8 ergs-i. 

The value of zint can be estimated then setting tc ^ ics, 
since the cloud should enter the jet before being signifi- 
cantly distorted by the impact of the latter. Once shocked, 
the cloud can suffer lateral expansion and conduction heat- 
ing, which can speed up fragmentation due to instabilities. 
In this work, we choose Zjnt fixing ics = '2tc: 



Zint ^ 5 X 10^^ ( "- \ 

VlO^cms ^/ 



10 



10, 



-1/2 



1044 erg s- 



1/2 



cm. 



(10) 



We note that the available power in the bow shock is Lbs ~ 
(i7c/ctj)Lj cx z~^. Therefore, the most luminous individual 
jet-cloud interaction will take place at z ~ Zint- 

The BLR size can be estimated through an empirical 
relation obtained from sources with a well stablished BLR, 
i.e. FR II radio galaxies. This relation is in general of the 
type i?bir oc ibij., where Lbir is the luminosity of the BLR 
and a ~ 0.5 - 0.7 (e.g. Kaspi et al. 2005, 2007; Peterson 
et al. 2005; Bentz et al. 2006). In this paper we use the 
following relations: 



i?blr -- 6 X 10 



16 



L 



blr 



1044 erg s-1 



and 

-Rblr 



2.5 X 10^^ 



L 



blr 



1044 erg s-1 



0.55 



(12) 



from Kaspi et al. (2005, 2007). 

In Fig. [3] we show the relation of zint and i?bir with Lj, 
assuming that Lbu- is a 10% of the disc luminosity, which is 
taken here equal to Lj. As seen in the figure, for reasonable 
parameters, the condition Zint < Rhh is fulfilled for a wide 
range of Lj. Figure |3] also shows the relation between zq and 
Mbh, which shows that for A/bh ^ 10^ Mq the jet could be 



Fig. 3. The interaction height, z-mt (red solid line), and the 
size of the BLR, i?bir (green dashed line), for different values 
of Lj (bottom horizontal axis). We have derived i?bir(ij) 
fixing Lbir = 0.1 Lj and plotted i?bir using Eqs. pT]) and 
(IT^ . In the same plot, the height of the jet base, zq (blue 
dotted line), is plotted as a function of Mbh (top horizontal 
axis). 

even not (fully) formed at the BLR scales at the lowest 
Lj -values. 

3. Non-thermal particles and their emission 

In the bow and cloud shocks, particles can be accelerated 
through diffusive shock acceleration (Bell 1978). However, 
the bow shock should be more efficient accelerating par- 
ticles than the shock in the cloud because Wbs 3> Wcs- In 
addition, the cloud shock luminosity is smaller than the 
bow-shock luminosity by ^ l/(2x"'^^). For these reasons, 
we focus here on the particle acceleration in the bow shock. 
In this section, we briefly describe the injection and evo- 
lution of particles, and their emission, remarking those as- 
pects that are specific to AGNs. The details of the emitting 
processes considered here (synchrotron, IC and pp interac- 
tions) can be found in Araudo et al. (2009) and references 
therein. 

First, one can estimate the non-thermal luminosity, Lnt, 
injected at zi^t in. the bow shock in the form of relativistic 
electrons or protons: 



(11) Lnt = VntL 



bs 



4xl0»(|li) 



Rr 



104-^cm 



(13) 



Li 



1044 erg s-i 



ergs 



Then, the accelerator/emitter magnetic field in the bow- 
shock RF (B) can be determined relating Ub = VBUnt, 
where Ub = B'^/Stt and Unt = int/(o'cc) are the magnetic 
and the non-thermal energy densities, respectively. For lep- 
tonic emission and to avoid supression of the IC channel, 
high gamma-ray outputs require jye well below 1. In this 
context, B can be parametrized as follows: 



B w 10 



VO.oJ V109cms-i/ V10i0cm-3j 



G. 



(14) 
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Regarding the acceleration mechanism, since the bow shock 
is relativistic and the treatment of such a shocks is com- 
plex (see Achterberg et al. 2001), we adopt the following 
prescription for the acceleration rate: 



E, 



O.lqBc, 



(15) 



similar to that in the relativistic termination shock of the 
Crab pulsar wind (de Jager et al. 1996). 

Particles suffer different losses that balance the energy 
gain from acceleration. The electron loss mechanisms are es- 
cape downstream, relativistic Bremsstrahlung, synchrotron 
emission, and external Compton (EC) and SSC. We note 
that B, Lnt and the accelerator/emitter size at z-mt are con- 
stant for different Lj and fixed 77B and ?7nt, and only Lbir 
and Ld are expected to change with Lj. Therefore, as long 
as the external photon fields are negligible, the maximum 
electron energy at Zint does not change for different jet pow- 
ers. 

In Fig. m the leptonic cooling timescales are plotted 
together with the escape time and the acceleration time for 
a bow shock located at Zint- A value for 773 equal to 0.01 has 
been adopted. The SSC cooling timescale is plotted for the 
steady state. The escape time downstream the relativistic 
bow shock is taken as: 



to. 



3 Re 



= lO-' 



Re 



10 



13, 



S . 



(16) 



Synchrotron and EC/SSC are the dominant processes in 
the high-energy part of the electron population, relativis- 
tic bremsstrahlung is negligible for any energy, and elec- 
tron escape is relevant in the low-energy part. This yields 
a break in the electron energy distribution at the energy 
at which the synchrotron/IC time and the escape time are 
equal. The Thomson to KN transition is clearly seen in the 
EC cooling curves, but is much smoother in the SSC case. 
The maximun electrons energy are around several TeV. 
Given the similar cooling timescale for electrons via rel- 
ativistic bremsstrahlung and protons through pp collisions 
(^broms/pp ~ 10^^/n s, whcrc n is the target density), pro- 
tons will not cool efficiently in the bow shock. Photomeson 
production can also be discarded as a relevant proton cool- 
ing mechanism in the bow shock due to the relatively low 
achievable proton energies and photon densities. The maxi- 
mum proton energy of protons is constrained making equal 
the acceleration time and the time needed to diffuse out of 
the bow shock. Assuming Bohm diffusion, idiff — 3 Z^/2 rg c 
(where r^ is the particle gyroradius), maximum proton en- 
ergy is: 



77" max 



O.lqBRc = 5x10^ 



B 
lOG 



Re 



101 



TeV. (17) 



Electrons are injected in the bow shock region following 
a power-law in energy of index 2.2 (Achterberg et al. 2001) 
with an exponential cutoff at the maximum electron energy. 
The injection luminosity is L^t- To first order, the electron 
evolution can be computed assuming homogeneous condi- 
tions, following therefore a one-zone approximation with 
all the mentioned cooling and escape processes. The for- 
mulae for all the relevant radiative mechanisms, as well as 
the solved electron evolution differential equation, can be 
found in Araudo et al. (2009). In some cases, SSC is the 
dominant cooling channel at high energies. In that case, 
the calculations have to be done numerically splitting the 




• ^1'- 


= 10" 


ergs" 


K,. 


= 10^ 


6 -1 
ergs 


■ L^ = 


= 10^' 


erg s"' 


K. 


= 10* 


ergs 



Synchrotron 



9 10 11 

log(E^ [eV]) 



Fig. 4. Acceleration gain (orange dot-dashed line), es- 
cape (blue dashed line) and cooling lepton timescales are 
plotted. SSC (green two-dot-dashed line) is plotted when 
the steady state is reached and EC for the both BLR 
(turquoise square line) and disc (maroon dotted line) pho- 
ton fields are shown for the conditions of faint (BLR: 10*^; 
disc: lO''^ erg s~^) and bright sources (BLR: 10"'^; disc: 
lO'*^ erg s~i). Synchrotron (red solid line) and relativistic 
Bremsstrahlung (violet dot-two-dashed line) are also plot- 
ted. We have fixed here t/b = 0.01. For protons, the pp cool- 
ing timescale is not shown for clarity, although it would be 
similar to that of relativistic bremsstrahlung. 



evolution of the electron population in different time steps. 
In each step, the radiation field is updated accounting for 
the synchrotron emission produced in the previous step un- 
til the steady state is reached. The duration of each step 
should be shorter for the earlier phases of the evolutionary 
state to account properly for the rise of the synchrotron 
energy density in the emitter. 

Once the steady-state electron distribution in the bow 
shock is computed, the spectral energy distribution (SED) 
of the non-thermal radiation can be calculated. The syn- 
chrotron self-absorption effect has to be taken into account, 
which will affect the low energy band of the synchrotron 
emission. At gamma-ray energies, photon-photon absorp- 
tion due to the disc and the BLR radiation are to be consid- 
ered, the internal absorption due to synchrotron radiation 
being negligible. Given the typical BLR and disc photon 
energies, ~ 10 eV and ~ 1 keV, respectively, gamma rays 
beyond 1 GeV and 100 GeV can be strongly affected by 
photon-photon absorption. On the other hand, for most of 
cases photons between 100 MeV and 1 GeV energies will 
escape the dense disc photon field. 

Although proton cooling is negligible in the bow-shock 
region, it may be significant in the cloud. Protons can pen- 
etrate into the cloud if tcsc > tdiff, which yields a minimum 
energy to reach the cloud of Ep ~ OAE^^^. These pro- 
tons will radiate in the form of gamma rays only a fraction 
- 3 X 10-1 (i?c/10^3cm)(ipp/10^s)^i of their energy, wich 
makes the process rather unefRcient. The reason is that 
these protons cannot be efficiently confined and cross the 
cloud at a velocity ^ c. For further details of the proton 
energy distribution in the cloud, see Araudo et al. (2009). 
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4. Many clouds inside the jet 

Clouds fill the BLR, and many of them can simultaneously 
be inside the jet at different z, each of them producing non- 
thermal radiation. Therefore, the total luminosity can be 
much larger than the one produced by just one interaction, 
which is ~ Lnt- The number of clouds within the jets, at 
z < i?bir, can be computed from the jet {Vj) and cloud (14) 
volumes, resulting in: 



^-2/| 



1044 



ergs' 



Rr 



10 



13, 



(18) 



where the factor 2 accounts for the two jets and / ^ 10^^ is 
the filling factor of clouds in the whole BLR (Dietrich et al. 
1999). Actually, N^ is correct if one neglects that the cloud 
disrupts and fragments, and eventually dilutes inside the 
jet. For instance, Klein et al. (1994) estimated a shocked 
cloud lifetime in several ics, and Shin et al. (2008) found 
that even a weak magnetic field in the cloud can signifi- 
cantly increase its lifetime. Finally, even under cloud frag- 
mentation, strong bow shocks can form around the cloud 
fragments before these have accelerated close to Vj . All this 
makes the real number of interacting clouds inside the jet 
hard to estimate, but it should be between (tcs/ij) N^ and 

The presence of many clouds inside the jet, not only at 
Zint but also at higher z, implies that the total non-thermal 
luminosity available in the BLR-jet intersection region is: 

az 



2 X 10« C^) 
\0.lJ 



Re 



10 



13, 



1044 erg s- 



1.7 



(19) 



where dN^/dz is the number of clouds located in a jet vol- 
ume dVj = 7r(0.1z)2dz. In both Eqs. [IHl and [H ibir has 
been fixed to 0.1 Lj, approximately as in FR II galaxies, 
and Rhh has been derived using Eq. ([TT|) . 

In Fig. [SI we show estimates for the gamma-ray lumi- 
nosity when many clouds interact simultanously with the 
jet. For this, we have followed a simple approach assuming 
that most of the non-thermal luminosity goes to gamma 
rays. This will be the case as long as the escape and syn- 
chrotron cooling time are longer than the IC cooling time 
(EC-I-SSC) for the highest electron energies. Given the lit- 
tle information for the BLR in the case of FR I sources, we 
do not specifically consider these sources here. 

In the two next sections, we present more detailed cal- 
culations applying the model presented in Sect. [3] to two 
characteristic sources, Cen A (FR I, one interaction) and 
3C 273 (FR II, many interactions). 



5. Application to FR I galaxies: Cen A 

Cen A is the closest AGN, with a distance d « 3.7 Mpc 
(Israel 1998). It has been classified as an FR I radio galaxy 
and as a Seyfert 2 optical object. The mass of the black 
hole is « 6 X 10^ M© (Marconi et al. 2000). The angle be- 
tween the jets and the line of sight is large, > 50° (Tingay 
et al. 1998), thus the jet radiation towards the observer 
should not suffer strong beaming. The jets of Cen A are 
disrupted at kpc scales, forming two giant radio lobes that 
extend ^ 10° in the southern sky. At optical wavelenghts. 
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Fig. 5. Upper limits on the gamma-ray luminosity pro- 
duced by N\ clouds inside the jet as a function of ij in 
FR II sources. Two cases are plotted, one assuming that 
clouds cross the jet without disruption (green solid lines), 
and one in which the clouds are destroyed in a time as short 
as tcs (green dashed lines). The thick (solid and dashed) 
and thin (solid and dashed) lines correspond to i?bir oc i^j^ 
(Kaspi et al. 2007) and i?bir oc Lgi^^ (Xaspi et al. 2005), 
respectively. In addition, the sensitivity levels of Fermi in 
the range 0.1-1 GeV (maroon dotted lines) are plotted for 
three different distances d = 10, 100 and 1000 Mpc. 



the nuclear region of Cen A is obscured by a dense re- 
gion of gas and dust, probably as a consequence of a recent 
merger (Thomson 1992, Mirabel et al. 1999). At higher 
energies, Chandra and XMM-Newton detected continuum 
X-ray emission coming from the nuclear region with a lu- 
minosity ~ 5 X 104^ erg s~^ between 2-7 keV (Evans et 
al. 2004). These X-rays could be produced by the accretion 
flow and the inner jet, although their origin is still unclear. 
In the GeV range, Cen A was detected above 200 MeV by 
Fermi, with a bolometric luminosity of « 4 x 104° erg s~^ 
(Abdo et al. 2009), and above - 200 GeV by HESS, with 
a bolometric luminosity of sa 3 x 10'^^ erg s~^ (Aharonian 
et al. 2009). In both cases, this HE emission is associated 
with the nuclear region. Cen A has been proposed to be a 
source of ultra HE cosmic rays (Romero et al. 1996). 

A BLR has not been detected so far in Cen A 
(Alexander et al. 1999), although this could be a conse- 
quence of the optical obscuration produced by the dust 
lane. One can still assume that clouds surround the SMBH 
in the nuclear region (Wang et al. 1986, Risaliti et al. 2002) 
but, as a consequence of the low luminosities of the ac- 
cretion disc, it is not expected that the photoionization of 
these clouds will be efficient enough to produce lines. Since 
no emission from these clouds is assumed, we only consider 
the EC scattering with photons from the accretion flow. 

We adopt here a jet power for Cen A of Lj = 
1044 erg s~^. From this value, and the values given in Sect. [2] 
for the remaining parameteres of the jet and the cloud, Zjnt 
results in Ri 5 X 10^^ cm. At this jet height, the emission 
produced by the interaction between one cloud and the jet 
is calculated assuming a t/b = 0.01, and the corresponding 
SED is presented in Fig.jB] As mentined in Sect. [31 the low- 
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Fig. 6. Computed SED for one jet-cloud interaction at Zint 
in Can A. We show also the SEDs of the detected emission 
by Fermi and HESS, as well as the sensitivity curves of 
these instruments. 

Table 2. Adopted parameters for Cen A and 3C 273. 



Fig. 7. Computed SED for one jet-cloud interaction at zint 
in 3C 273. The emission in the 0.1-1 GeV range from many 
clouds inside the jet is also shown, together with the sensi- 
tivity level of Fermi and the observed SED above 200 MeV. 



Gen A 



3G 273 



Distance [Mpc] 
Black hole mass [Mq] 
Inclination angle [°] 
Jet luminosity [erg s~^] 
Disc luminosity [erg s~^] 
Disc photon energy* [eV] 
BLR luminosity Lbir [erg s^^] 



3.7 

6 X 10^ 

> 50 

^q44 
5 X W^ 

~ 5 X 10^ 



6.7 X 10^ 
7x 10^ 

~ 15 
4 X lO"*^ 
2 X 10""^ 

54 
4 X lO"*^ 



of the thermal component. 



energy band of the synchrotron spectrum is self-absorbed at 
energies below ~ 10~ eV. At gamma-ray energies, photon- 
photon absorption is negligible due to the weak ambient 
photon fields (e.g. Rieger & Aharonian 2009, Araudo et al. 
2009, 2010a, b). At high energies, SSC dominates the radia- 
tive output, with the computed luminosity above 100 MeV 
being ~ 2 x lO^^ erg s-\ and above 100 GeV about 10 
times less. These luminosities are below the sensitivity of 
Fermi and HESS and one order of magnitude smaller than 
the observed ones. Note however that Lnt c>c R^, and for 
slightly bigger clouds, int may grow up to detectable lev- 
els. The penetration of a big clump in the base of the jet of 
Cen A would lead to a flare with a duration of about one 
day. 



6. Application to FR II galaxies: 3C 273 (ofF-axis) 

3C 273 is a powerfuU radio-loud AGN at a distance of 
d = 6.7 X 10^ Mpc (Courvoisier 1998) with a SMBH mass 
Mbh -- 7 X lO^M© (Paltani & Tiirler 2005). The angle of 
the jet with the line of sight is small, ~ 6°, which implies 
the blazar nature of 3G 273 (Jolley et al. 2009). The whole 
spectrum of this source shows variability (e.g. Plan et al. 
1999) from years (radio) to few hours (gamma rays). At 
high energies, 3C 273 was the first blazar AGN detected 
in the MeV band by the GOS-B satellite and, later on, by 



EGRET (Hartman et al. 1999). Recently, this source was 
also detected at GeV energies by Fermi and A GILE, but it 
has not been detected yet in the TeV range. Given the jet 
luminosity of 3C 273, Lj w 4 x 10^^ erg s""'^ (Kataoka et al. 
2002), Zint results in « 3 x 10^'' cm. The BLR luminosity of 
this source is « 4 x lO''^ erg s~^ (Cao & Jiang 1999), and its 
size 7 X 10^^ cm (Ghissellini et al. 2010), which implies that 
jet-cloud interactions can take place. The disc luminosity 
is high, « 2 X 10^^ erg s~^, with typical photon energies 
« 54 eV (Grandi & Palumbo 2004). 

The non-thermal SED of the radiation generated by 
jet-cloud interactions in 3C 273 is shown in Fig. [71 At 
Zint, the most important radiative processes are synchro- 
ton and SSC. The bolometric luminosities by these pro- 
cesses in one interaction at Zjnt are 6 x 10'^* erg s~^ and 
2 X lO'^® erg s~^, respectively. Given the presence of the 
strong radiation fields from the disc and the BLR, the emis- 
sion above ~ 10 GeV is absorbed through photon-photon 
absorption, and the maximum of the emission is around 
0.1-1 GeV. Given the estimated number of clouds in the 
BLR of 3C 273, - 10^ (Dietrich et al. 1999), and the size is 
R^jIj. w 7 X 10^^ cm, the filling factor results in / ^ 3 x 10"^. 
With this value of /, the number of clouds in the two 
jets results in ~ 2 x 10'^ and 5 x 10^ for both the mini- 
mum and the maximum values (see Sect.0]). Considering 
the most optimistic case, the SSC luminosity would reach 
2 X 10^''' erg s~^. This value is well below the observed lu- 
minosity by Fermi in the GeV range, ~ 3 x 10'*^ erg s^^ in 
the steady state and ~ 1.7 x 10^^ erg s~^ in flare (Soldi et 
al. 2009). However, the detected emission is very likely of 
beamed origin and should mask any unbeamed radiation. 
However, powerful non-blazar AGNs (FR II galaxies) do 
not present this beamed component, which makes possible 
the detection of GeV emission from jet-cloud interactions in 
these sources. In this case, given that many BLR clouds can 
interact with the jet simultaneously, the radiation should 
be steady. 



Araudo, Bosch-Ramon & Romero: Gamma rays from jet-cloud interactions in AGNs 



7. Summary and discussion 

In this work, the interaction of clouds with the base of jets 
in AGNs is studied. Considering reasonable cloud and jet 
parameters, we estimate the relevant dynamical timescales 
of these interactions, concluding that clouds can enter into 
the jet only above a certain height, ~ z-mt- Below z-mt, the 
jet is too compact and its ram/magnetic pressure will de- 
stroy the cloud before fully penetrating into the jet. Once 
the cloud significantly interacts with the jet, strong shocks 
are generated and gamma rays can be produced with an ef- 
ficiency that depends strongly on the bow-shock magnetic 
field. 

Bow shock _B-values well below equipartition with non- 
thermal particles allow significant gamma-ray emission. For 
very high B- values (Poynting-flux dominated jets), the 
treatment performed here does not apply. In that case, Zint 
could be still defined adopting the jet magnetic pressure 
instead of the ram pressue. If a cloud entered in such a jet, 
particle acceleration in the bow shock could still occur due 
to, for instance, magnetic reconnection. The study of such 
a case would require a completely different approach than 
the one presented here. In general, for bow-shock magnetic 
fields above equipartition with non-thermal particles, the 
IC channel and gamma ray production will be supressed in 
favor of synchrotron emission. Unless magnetic dissipation 
reduced the magnetic field enough for IC to be dominant. 
Bow shock i3-values well below equipartition with non- 
thermal particles allow significant gamma-ray emission. We 
note that modeling of gamma rays from AGN jets uses to 
require relatively low magnetic fields (e.g. Ghisellini et al. 
2010). Therefore, it could be that, even if the jet magnetic 
field were high at z-mt, it could become small enough farther 
up due to bulk acceleration (e.g. Komissarov et al. 2007) or 
some other form of magnetic dissipation. 

For very nearby sources, like Cen A, the interaction of 
big clouds with jets may be detectable as a flaring event, al- 
though the number of these big clouds and thereby the duty 
cicle of the flares are difficult to estimate. Given the weak 
external photon flelds in these sources, VHE photons can 
escape without suffering significant absorption. Therefore, 
jet-cloud interactions in nearby FR I may be detectable in 
both the HE and the VHE range as fiares with timescales 
of about one day. Studying such a radiation would provide 
information on the environmental conditions and the base 
of the jet in these sources. 

In FR II sources, many BLR clouds could interact si- 
multaneously with the jet. The number of clouds depends 
strongly on the cloud lifetime inside the jet, which could be 
of the order of several tcs- Nevertheless, it is worthy noting 
that after cloud fragmentation many bow shocks can still 
form and efficiently accelerate particles if these fragments 
are slower than the jet. Since FR II sources are expected 
to present high accretion rates, radiation above 1 GeV pro- 
duced in the jet base can be strongly attenuated due to 
the dense disc and the BLR photon fields, although gamma 
rays below 1 GeV should not be significantly affected. Since 
jet-cloud emission should be rather isotropic, it would be 
masked by jet beamed emission in blazar sources, although 
since powerful/nearby FR II jets do not present significant 
beaming, these objects could indeed show gamma rays from 
jet-cloud interactions. In the context of AGN unification 
(Urry & Padovani 1995), the number of non-blazar (radio- 
loud) AGNs should be much larger than that of blazars 



with the same Lj. As shown in Fig. [SJ close and powerful 
sources could be detectable by deep enough observations of 
Fermi. After few-years exposure a significant signal from 
these objects could arise, their detection providing strong 
evidence that jets are already strongly matter dominated 
at the bow shock regions, as well as physical information 
on the BLR and jet base region. 
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